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Abstract. A fraction of neutrino emission from GRB accretion disks annihilates above the disk, creating e+- plasma 
that can drive GRB explosions. We calculate the efficiency of this annihilation using the recent detailed model of 
hyper-accretion disks around Kerr black holes. Our calculation is fully relativistic and based on a geodesic- tracing 
method. We find that the efficiency is a well-defined function of (1) accretion rate and (2) spin of the black hole. 
It is practically independent of the details of neutrino transport in the opaque zone of the disk. The results help 
identify the accretion disks whose neutrino emission can power GRBs. 
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The relativistic model of GRB accretion disks w^as completed recently fT\. It describes the disk matter dov\^n 
to the last stable orbit r^^ and gives energy fluxes carried away by v and v at all radii r. In the present w^ork, 
we trace the trajectories of emitted neutrinos, and calculate the rate of vv annihilation around the disk. This 
process deposits plasma and can play a key role in the formation of GRB jets. 

Neutrino annhilation was previously calculated in a number of w^orks (e.g. O, ||3l, ||4|) Our w^ork has 
two motivations: (1) A relativistic calculation has never been done for a realistic accretion disk around a 
spinning black hole. Previous works either used a toy model for neutrino source (e.g. an isothermal disk, 
HI) or replaced neutrino trajectories by straight lines |[2|. (2) The efficiency of vv annihilation depends 
strongly on the accretion rate M and the spin parameter a of the black hole. It is desirable to know this 
dependence and identify the range of M and a where vv annihilation can provide the observed energy of 
GRB explosions. 



Fortunately, our final result depends only on the energy fluxes Fy and Fy from the disk surface, which are 
independent of the neutrino-transport details and are already calculated in The rate of vv annihilation is 
insensitive to the exact shapes of v and v spectra. This fact can be demonstrated using two extreme models 
A and B: 

Model A: Neutrinos v and v are emitted with the same spectrum as found inside the disk (same 
temperature T and chemical potential jiy). The spectrum is normalized so that the emerging emission carries 
away the known energy fluxes Fy and Fy. 

Model B: Neutrinos are emitted with the effective surface temperature T^ff defined by {7/^)aT^^ = Fy-\-Fy 
(a is the Stefan-Boltzmann constant). The two models give practically the same vv annihilation rate 
(Figure 2). Our results confirm the analytical argument that it is sufficient to know T^ff to calculate neutrino 
annihilation rate [5J. When the disk is efficiently cooled (neutrino losses almost balance viscous heating), 
Teff is given by the standard thin-disk model [6 J: Tgff = rj^^ndard j^^^^ model applies to GRB disks in a broad 
range of accretion rates Mign <M < Mtrap E/ where 
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NEUTRINO EMISSION FROM THE DISK 
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FIGURE 1. Spatial distribution of the energy deposition rate by vv annihilation. This example assumes the accretion 
rate M = \Mqs~^, the disk viscosity parameter a = 0.1, and the black hole mass M = 3Mq; Model A is assumed for 
the neutrino spectrum (see the text). Left panel: a = 0. Right panel: a = 0.95 {a= 1 corresponds to the maximally rotating 
black hole). Note that the energy deposition rate is much higher in the case of a = 0.95. The arrows show the specific 
momentum of the plasma injected by vv annihilation. The white curve is where the radial component of the injected 
momentum changes sign. This boundary gives an idea of the region where the injected plasma is lost into the black 
hole. 



Here a ~ 0. 1 is the standard viscosity parameter of the accretion disk, and the factors K depend on the black 
hole spin a; e.g. for a = 0.95 they are ^ign = 0.021 M© s~^ and ^trap = 1.8 M© s~^ T^ff of neutrino emission 
can be approximately described as 

r M<Mign 
Teff (M, r) ^ 7;r^^^d(Mign, r) x \ (M/Mign)^/' <M< Mrap (2) 

[ (Mtrap/Mign)'/' M > Mtrap 

Teff is suppressed for M < Mign because the disk is not hot enough to ignite neutrino-emitting reactions. 
Teff suddenly increases when M = Mign and growls as M^/^ until M = Mtrap- At higher M, neutrinos become 
trapped in the disk (advected into the black hole) and T^ff saturates. Equation (|2| defines our Model C, 
w^hich reproduces surprisingly well the more detailed numerical results (Figure 2). This model allows us to 
obtain an explicit approximate formula for the annihilation rate (eq.lslbelow). 



RESULTS 

Figure 1 shows two examples of the spatial distribution of the energy deposition rate by vv annihilation. 
Integration of this distribution over volume outside the black hole gives the total energy deposition rate 
Eyy. We performed this calculation for various M (Fig. 2). Our results show that Eyy is well approximated 
by a simple formula, 

r M<Mign ] 

£vv ^ 9 X 10^1 X't-^ X <^ m9/4 Mign <M< Mtrap \ ^ (3) 
I ^top ^>Mrap J 

where m = M/MqS~^ and Xms = rnis(<3)(2GM/c^)~^ Derivation of the scaling of Eyy with m and x^s is given 
in |5|. The dependence of Eyy on the black hole spin is huge: x^f^ varies by a factor of 170 for < a < 0.95. 
Note that a (viscosity parameter of the disk) enters the result only through Mign and Mtrap (eq.[l|. 



2 3 4 5 6 

FIGURE 2. Left panel: Total energy deposition rate due to vv annihilation outside the black-hole horizon, E^^, as 
a function of the disk accretion rate, M. The two characteristic accretion rates Mig^ and Mtrap depend on the viscosity 
parameter a (see eq.[l|; a = 0. 1 is chosen in this figure. The black hole is assumed to have mass M = 3 M©. Eyy strongly 
depends on the spin parameter of the black hole; the results are shown for two cases: a = (triangles) and a = 0.95 
(squares). The uncertainty in the vertical structure of the accretion disk leads to a small uncertainty in Eyy as illustrated 
by two extreme models: Model A (open symbols) and Model B (filled symbols), see the text for details. The results of 
both models are well approximated by Model C (eq.jsj, which is shown in the figure by line; the line is dotted at low M 
where the disk is transparent to neutrinos. Right panel: Dependence of the energy deposition rate Eyy on the black hole 
spin for a fixed M = \Mq/s. Instead of using the spin parameter a directly, it is more instructive to plot Eyy versus the 
radius of the last (marginally stable) orbit rms- Then one can see the power-law dependence of Eyy on rms^ Eyy oc r~g 
The radius rms is a function of a, for non-rotating black holes rms = 6GM/c^ and for maximally rotating black holes 
rms = GM/c^. 
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The efficiency of vv annihilation can be defined as e = Eyy /L where L is the total neutrino luminosity of 
the disk. For example a = 0.95 (which corresponds to x^s ^ 1) gives L ^ O.lMc^ and 

e~0.05f ^j^— , Mign<M<Mttap. (4) 

The observed GRB luminosity Lobs can be supplied by vv annihilation around a rapidly spinning black 
hole (a = 0.95) if M > O.38(M0/^) (Lobs/10^' erg/s), which is in the range of plausible accretion rates in GRB 
central engines. 

Finally, note that Eyy is defined as the total energy deposition rate outside the event horizon. A fraction 
of the created plasma falls into the black hole and not contribute to the observed explosion (Figure 1). 
The corresponding refinement of £ depends on the plasma dynamics outside the disk, which is affected by 
magnetic fields and is hard to calculate without additional assumptions. 
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